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Energy fluctuations in a biharmonically driven nonlinear system 

Navinder Singh, Sourabh Lahiri and A. M. JayannavaiQ 
Institute of Physics, Sachivalaya Marg, Bhubaneswar-75 1 005 , India 

We study the fluctuations of work done and dissipated heat of a Brownian particle in a symmetric 
double well system. The system is driven by two periodic input signals that rock the potential 
simultaneously. Confinement in one preferred well can be achieved by modulating the relative phase 
between the drives. We show that in the presence of pumping the stochastic resonance signal is 
enhanced when analyzed in terms of the average work done on the system per cycle. This is in 
contrast to the case when pumping is achieved by applying an external static bias, which degrades 
resonance. We analyze the nature of work and heat fluctuations and show that the steady state 
fluctuation theorem holds in this system. 

PACS numbers:05.40.-a; 05.40. Jc; 05.60. Cd; 05.40. Ca 
Key Words: Stochastic Resonance, Fluctuation Theorem 

I. INTRODUCTION 

A number of important results have been obtained over the last two decades on the statistical 
properties of fluctuations in physical quantities in non-equilibrium processes. These are referred 
to as Fluctuation Theorems (FTs) 0, 0, |, fl, H, 0, 0, 0, E EQ tH, EE El, EE EE E3- They 
allow extension of thermodynamic concepts to small systems [if. The fluctuation theorems reveal 
rigorous relations for the properties of distribution functions of physical variables such as work done, 
dissipated heat or entropy production for systems driven far away from equilibrium, independent 
of the nature of driving. They are not restricted to the linear response regime, thus allowing us 
to obtain results generalizing Onsager reciprocity relations to the nonlinear response coefficients 
in nonequilibrium state. From these theorems, corollaries such as the statistical derivation of the 
second law, can be established. There are different fluctuation theorems, depending on the physical 
quantities they relate to and on the state of the system they refer to. These theorems are useful 
to probe nonequilibrium states in nanophysics and biology. Hence they are anticipated to play an 
important role in the design of nanodevices and engines (molecular motors) at nanoscales. The 
distributions of heat and work in relation to FTs have been experimentally studied for few Brownian 
systems 0, El, M, 



In recent theoretical [24j . |25j and experimental [2 



work done on the system have been explored in a system exhibiting stochastic resonance [24], [25 



studies, the distributions of dissipated heat and 

The 



steady state fluctuation theorem (SSFT) (as elaborated upon in section IIIIFj) holds in this system. 
Exploring the FTs in nonlinear systems by changing the symmetry of the driving force cycle has been 
suggested in Ref . [26[ . To this end, we study the dynamics of a particle in a symmetric double well 
potential which is in contact with a thermal bath at temperature T. This system exhibits stochastic 
resonance (SR) under subthreshold external ac drive [27|. The fundamental periodic component of 
the system response (i.e., the amplitude of the feeble input at the same frequency) can be amplified 
by the assistance of noise. It is reflected as a peak in the output signal-to-noise ratio as a function 
of noise strength. This peak occurs when the noise induced switching rate in the system matches 
the forcing frequency. This optimization or synchronization condition is achieved by tuning the 



noise intensity. This phenomenon is known as stochastic resonance [271 . [28|. Here noise plays a 
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constructive role as opposed to our conventional wisdom that the presence of noise degrades the 
signal. This is due to a cooperative interplay between the system nonlinearity and input signal in 
addition to the noise. Because of its generic nature, this phenomenon boasts applications in almost all 
areas of natural science [28|. To characterize this resonance phenomenon, several different quantifiers 
have been introduced in the literature 28, m, m, eh, m, l33l . 0, [35|]. One of the quantifiers, namely 
the input energy of the system or the work done on the system per cycle is known to characterize 



SR as a bona fide resonance [2J, [25J, [3J, 135J. In this case, the resonance can be shown to occur both 
as a function of noise strength and driving frequency. 
It is known that static asymmetry in the bistable potential weakens the magnitude of the SR effect 

MM 



Static tilt in the potential makes one potential well more stable than the other leading to 
more particle localization or pumping in one well (lower well) compared to the other. Moreover, due 
to asymmetry in the potential, escape rate of a particle from higher to lower well will be different 
from lower to higher well. These two different rates make synchronization difficult between the 
signal and the dynamics of the particle hopping, since the driving frequency cannot match both 
these hopping rates simultaneously. 

In the present work, we study the SR for a particle in a symmetric double well potential, driven 
simultaneously by two periodic signals of frequencies uj and 2u with a relative phase difference 
between them. Such a force averaged over a period does not lead to a net bias and yet particle is 
jreferentially p umped i nto one well depending on phase difference <fi and other ph ysic al parameters 

Due 




0, H H El, H EH, El . This phenomenon is known as harmonic mixing 

to this statistical confinement of the particle, similar to the case of static tilt [28|, |43|J, we expect to 
observe a reduced SR signal in this system. However, contrary to this expectation, we show that the 
resonance signal is enhanced in the presence of the biharmonic drive at frequency 2a; when analyzed 
in terms of the input energy (or the work done) as a quantifier of SR. Using stochastic energetics 



13l . |4J, |45j we also study the nature of fluctuations in the work done, dissipated heat and internal 



energy across SR. In some range of parameters, nature of hysteresis loops is analyzed. We show that 
the SSFT holds for work done over a long time interval. The modified SSFT for heat is also studied 
in this system. 



II. THE MODEL: BROWNIAN PARTICLE IN A ROCKED DOUBLE WELL POTENTIAL 

We consider the stochastic dynamics of a Brownian particle in a double- well potential V(x) = —\-\- 

rocked by a weak biharmonic (time-asymmetric) external field F(t) = A cos(tut) + B cos(2u>t + 0). 
The potential V(x) has two minima at x = ±1, separated by a central (at x = 0) potential barrier 
of height AV = 0.25. The overdamped Langevin dynamics is given by [461 ] . 

where U(x,t) = V(x) — xF(t), 7 is the friction coefficient, £(t) is the Gaussian white noise with 
the properties 

(tit)) = 0, 

(mm) = 2DS(t-f), (2) 

where D = jknT. The thermodynamic work done by an external drive over a period r £J (= — ) is 



riven by [13|, |44], |45 
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fto+r^ dUix.t) , 

W > = L dt 



/ x(t)[Ausmujt + 2Bu;sm(2ujt + (f))]dt. (3) 

Jt 



This work (or input energy) over a period equals the change in the internal energy AU P = U (a; (to + 
Tuj),to + rj) — U(x(to),t Q ) plus the heat dissipated over a period Q p , i.e., 

W p = AU P + Q p . (4) 
The above equation is the statement of the First law of thermodynamics and can readily be 



obtained using stochastic energetics [44J. Since x(t) is a stochastic variable, it follows from eq. (3) 
and eq.(4) that W p , AU P and Q p are random variables when evaluated over different periods and 
realizations of x(t). The above model is solved numerically by using Heun's method [47] (all the 
physical quantities are in dimensionless units). We have ignored the initial transient regime up to 
time to and evaluated W p , Q p , and AU P over many cycles (~ 10 5 ) of a single long trajectory of the 
particle. 

III. RESULTS AND DISCUSSIONS 
A. SR as a function of noise strength 

In figure [TJ we have plotted the average work done over a single period (W p ) in the time asymptotic 
regime as a function of noise strength D for different values of biharmonic drive strength B (for 
A=0.1). Phase difference is taken to be zero. Other parameters are mentioned in the figure captions. 



For the case of B = we have reproduced earlier results [24J, |25|, [34j, [35] . The average input energy 
(Wp)) shows a peak signifying stochastic resonance (SR) as discussed extensively in earlier literature 
24, 25, 34, 35[. The quantity (W p ) can also be identified as the average dissipated heat or hysteresis 
loss into the bath in a time periodic steady state. This follows from eq. (4) by noting that the 
internal energy being a state variable, AU P averaged over a period is identically equal to zero. As 
we increase B, (W p ) increases for fixed value of D. Moreover, for different values of B, the system 
exhibits SR as a function of noise strength. The system in a steady state absorbs energy from the 
external drive and the same is dissipated as heat into the surrounding medium. Hence it is expected 
that at the resonance the system will absorb maximum energy from the external drive. The input 
energy curves for higher values of B lie above those for the lower values of B. The graph for B = 



matches with the earlier known results [25], |34|. |35| . and the peak position shifts towards higher values 



of D. It is evident from the figure that in the presence of biharmonic drive enhancement of SR signal 
occurs even though there is more statistical confinement of the particle (as B increases) in one well 
as shown in figure 2. In this figure we have plotted average position ((x)) over period in the time 
asymptotic regime as a function of B for fixed D = 0.05. The value of (x) not being zero signifies 
selective pumping or localization of particle from one well to another in the presence of biharmonic 
drive. Correspondingly, the probability density distribution of the particle averaged over a period 
shows a marked asymmetry even though the potential V(x) is symmetric [38]. In the absence of 
second harmonic component i.e., B = 0, (x) = as expected. The pumping is very significant at 
low values of temperature. As we increase temperature, the effective pumping reduces. Around and 
beyond SR, pumping is quite small as shown in the inset of figure HJ 



4 



0.15 




FIG. 1: The input energy (W p ) as a function of D for different values of the strength of second harmonic (B). The parameters 
are: ui = 0.1, A = 0.1, and (f> — 0. 




FIG. 2: Particle mean position (a;) as a function of B for D = 0.05. In the inset we have plotted (x) as a function D. Other 
parameters are: B — 0.06, A — 0.1, and to — 0.1. 



Stochastic resonance being a synchronization phenomenon [30j, |33| it is expected that particle 
hopping dynamics between the wells get synchronized with the input signal. We find that the relative 



variance (RV) in physical quantities such as work 

2fJ. 



show minima at SR 1 24, 25 



(w P ) 



and heat 



(Q P ) 



also 



In figure [3] we have plotted relative variance (RV) as a function of D for various values of B. The 
parameters used are the same as in figure 1. For a given value of B the RV shows a minimum around 
the same value of D at which (W p ) exhibits a maximum. As the amplitude of the biharmonic drive 
B increases, RV curves shift downwards. Higher the value of B, the lower is the value of RV at 
the resonance. These results are consistent with figure 1. In the parameter regime that we have 
considered, the RV is larger than one, i.e., variance in work is large compared to the mean. Hence 
in this regime, one should analyze full probability distribution as opposed to moments to get better 
understanding of the phenomenon. 

Increasing the amplitude of biharmonic drive leads to more statistical confinement of parti- 
cles (figure [2]). This must be reflected in the nature of hysteresis loops [H, 40]. More the pumping, 
more is the asymmetry in the hysteresis loops, as can be seen in figure lU^a) and (b). In these figures 



0.1 0.2 0.3 0.4 

D 



FIG. 3: Relative variance (RV) of input energy versus D for different values of B. Other parameters are same as in figure 1. 




F(t) F(t) 

(a) (b) 

FIG. 4: (a) Hysteresis loops ((a;)) versus F for different values of B, and for D — 0.1. Other parameters are same in figure 1. 
(b) Hysteresis loops for different values of B at D — 0.05. 

hysteresis loops are plotted for different values of D and B. The pumping of the particles also gets 
reflected in the shifting of figures in the vertical upward direction(as (x) > 0). For the case when 
B — 0, there will be no pumping and as expected loop is symmetric. 

B. SR in the presence of static tilt 

Particle pumping in a preferential well can also be induced by applying a static tilt to the otherwise 
symmetric double well potential. For this we take potential to be V\(x) = —\ + \ — ex. Depending 
on the value of c, the potential V\{x) becomes asymmetric and obviously more pumping results in 
the lower potential well. When this system is driven by external AC force A cos out we show that SR 
signal weakens. Figure GD^a) shows the average input energy as a function of D for various values of 
c. From this, we notice that the input energy curves for higher value of c are below those with lower 
value of c (other parameters being fixed). As c increases SR peak becomes broadened and shifts 
towards higher values of D. We thus observe that in the presence of pumping induced by static tilt, 



0.2 
D 
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D 



(a) 



(b) 



FIG. 5: (a) Plots of (W p ) as a function of D for different values of the static tilt (c). (b) Corresponding plots of relative variance 
of W p as a function of D. Fixed parameters are mentioned on the graphs. 
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FIG. 6: Plots showing hysteresis loops for different values of static tilt (c). Fixed parameters are: D = 0.1, A = 0.1 and u — 0.1. 

SR weakens as mentioned in the introduction. This is also corroborated by the nature of relative 
work fluctuations as a function of D (figure E(b)). From this figure we note that as we increase c 
RV increases for a given value of D. The magnitude of the RV at the minimum becomes larger as 
we increase c. This implies degradation of SR signal in the presence of particle pumping induced 
by a static tilt. The pumping due to static tilt makes the hysteresis loops asymmetric (figure [6]). By 
increasing c more pumping is achieved and this is reflected in the vertical shift of hysteresis loops. 
Thus from the above figures and discussions, we conclude that in the presence of biharmonic drive, 
SR increases while in the presence of static tilt, SR weakens. 



C. SR as a function of driving frequency 



In figure [7](a), we have plotted average input energy as a function of u for various values of B. 
Once again we notice that SR signal even for this case is increased as we increase the biharmonic 
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(a) (b) 

FIG. 7: (a) The average input energy per period (W p )as a function of frequency uj. (b) Relative variance (RV) of input energy 
vs frequency uj for different values of the strength of the second harmonic, B. The parameters used are: u — 0.1, A — 0.1, (f> — 0. 



component B. Each curve exhibits a peak as a function of u, thus establishing SR as a bona fide 
resonance. The peak shifts to the lower values of to as we increase B. This is consistent with the fact 
that peaks in figure [1] shift towards larger values of D as we increase B. This is a requirement for the 
time scale matching between D and uj. Since increase in B slows down the effective time averaged 
hopping rates between the wells, hence higher D is required to achieve resonance. This lowering of 
effective escape rate at given D in turn implies decrease in the resonant frequency. The enhancement 
of SR signal in the presence of B can be inferred from figure [3(b) where we have plotted relative 
variance across the SR as a function of uj for various values of B. Lower values of relative variance 
across the SR for larger values of B are suggestive of the fact that SR is enhanced in the presence 
of biharmonic drive, consistent with the conclusions of figure [TJ 

D. Energy fluctuations over a single period 

Next, we analyze the nature of distribution functions of input energy P(W P ), dissipated heat P(Q P ) 
and internal energy P(AU p ) for different values of D. These distributions are plotted in figured] (a), 
(b), and (c) below resonance (D = 0.05), at resonance (D = 0.12), and above resonance (D = 0.3) 
respectively. The averaged internal energy (U) being a state function assumes the same value at 
the beginning and at the end of a period or periods in the time asymptotic regime. Hence total 
change in the internal energy (AU P ) average over a period is equal to zero and it is also expected 
that the distribution P( AU„ ) is symmetric as is evident from figure [8](a),(b), and (c). The nature of 
P(AU P ) is explained in [26| for a single harmonic drive. As opposed to AU P , distributions for W p 
and Q p are asymmetric. These distributions keep on changing in shape depending on the number of 
cycles over which they have been obtained which will be discussed later in connection with steady 
state fluctuation theorem (SSFT). Probability distributions for work and heat have finite weights 
for the negative values of their arguments. These negative values correspond to the trajectories 
where the particle moves against the perturbing AC field over a short time. For small values of 
D(D = 0.05), peak for W p or Q p near the origin corresponds mainly to intrawell dynamics of the 
particle and is mostly confined to a single well. The occasional excursion of the particle into the 
other well as a function of time is clearly reflected as a small hump at higher values of W p or Q p in the 
plot of P(W P ) and P(Q P ). As we increase D further (D = 0.12), interwell dynamics starts playing 
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FIG. 8: Plots (a), (b), and (c) show the distributions P(W P ), P(Q P ), and P(AU P ) for different values of D, below resonance 
D = 0.05, at resonance D — 0.12, and above resonance D — 0.3, respectively. Other fixed parameters are also shown on the 
graphs. 

a dominant role, and hence the distributions become broader. Work distribution exhibits three 
prominent peaks including one at the negative side. For larger values of D beyond SR point, shapes 
of P{W P ) and P{Q P ) tend closer to Gaussian distribution with increased variance/fluctuations. For 
such high temperatures, particle makes several random excursions between the two wells during a 
single time period of the external drive. It may be noted that the relative variances in W p and Q p 
are larger than 1. Also, fluctuations in heat are larger than those of work when averaged over a 
single period. 

E. Effect of phase difference on SR 

We now analyze the role of phase difference (0) between driving fields on pumping and energetics 
of the system. In figure [9] (a), we have plotted (W p ) as a function of noise strength D for various 
values of <fi. Other physical parameters are held fixed as mentioned in the figure captions. In figure 
M (b), we have plotted relative variance of W p as a function of D. It is interesting to note that (W p ) 
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(a) (b) 

FIG. 10: Average position (a;) as a function of phase <j> for two different values of D. In (a) D = 0.05, and in (b) D = 0.12. 
Other fixed parameters are: B = 0.06, A = 0.1, and to — 0.1 



is insensitive to 0, even though the relative variance depends on 0. This is a rather surprising result, 
given the fact that different values of phase <fi lead to different degrees of localization of the particle 
in one of the wells. We have characterized this dynamic localization of particles by average position 
(x) of the particle in the double well potential which in fact can be large depending on D and <fi. This 
is shown in figure [10] where we have plotted (x) as a function of <fi for two different values of noise 
strength D. One can readily see that (x) is periodic in as expected. The insensitivity of (W p ) on 
phase gets reflected in the hysteresis loop areas as shown in figures fTTl (a) (D = 0.1) and (b)(D = 0.05) 
for different values of and fixed value of B(B = 0.06). We notice that the areas of the hysteresis 
loops remain same for different 0. However, their shapes are asymmetric and qualitatively different 
for different 0(i.e., sensitive dependence on phase <p). Due to the different degree of localization or 
pumping, loops are shifted in (x) — F plane. The sensitivity of full probability distribution on the 
phase difference can be seen from figures [12j In these figures we have plotted P(W P ) and P{Q P ) 
for different values of <fr as indicated. Note that the distributions exhibit qualitative differences for 
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FIG. 11: (a) Hysteresis loops for different values of <j> at D — 0.1, (b) Hysteresis loops for different cj> a t D — 0.05, with other 
parameters A = 0.1, B — 0.06, and w = 0.1 




(a) (b) 



FIG. 12: Figures (a) and (b) show the distributions P(W V ) and P(Q P ) respectively for three different values of <f> as mentioned 
on the graphs. Here D — 0.12 and B — 0.06. 



different <\>. We have also verified separately that for different values of rocking amplitudes, as long 
as we are in subthreshold regime, average input energy is not very sensitive to cf) as opposed to 
full probability distribution and hysteresis loops. By tuning <f), one can achieve different degrees of 
particle confinement and can control the fluctuations in heat and work. 



F. Energy fluctuations and SSFT 



Finally we discuss the validity of SSFT in the present case of nonequilibrium time periodic 
steady state. SSFT implies the probability distribution of physical quantity A to satisfy relation 
p(A)/p(— A) = e^ A , where (5 is the inverse temperature of the bath and A is the work done on the 
system or the heat released to the bath over a long time of observation. For linear driven systems 
SSFT for work is indeed satisfied even for work done over a single period in the time asymptotic 



regime \4m . However, for nonlinear systems it has been observed experimentally and theoretically 
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FIG. 13: (a) The evolution of P(W P ) over different periods. In the inset P(Wio P ) is plotted together with its Gaussian fit 
G(W). (b) The plot of symmetry function (In p^L^^ ) versus j3W for various values of periods. The parameters used are: 
D — 0.16, uj — 0.1, B — 0.06, A = 0.1, = 0. The solid line is the best fit for symmetry function calculated for 10 cycles. 




(a) (b) 



FIG. 14: (a) The evolution of P(Q P ) over different periods. In the inset P(Qio p ) is plotted together with its Gaussian fit 
G(Q). (b) The plot of symmetry function (In p^Q\ ) versus f3Q for various values of periods. The solid line is the best fit for 
symmetry function calculated for 10 cycles. The parameters used are same as in figure 13. 



that SSFT is satisfied if one considers work done over a large numbers of cycles 24, 25, 26[. In regard 
to heat, SSFT is known to be valid for Q < (Q). Since (Q) increases with the number of periods 
or measured time interval in the limit of large n (n — > oo), (Q) — > oo and hence the conventional 
SSFT is valid over an entire range of Q j23[. It may be noted that there exists an alternative rela- 
tion for heat fluctuation, namely extended heat fluctuation theorem [l5|, 16[ . In figure [13] we have 
plotted probability distribution P(W np )of work W np integrated over different number (n) of periods. 
P(W np ) for a single period exhibits double peak structure. As we increase the number of periods the 
probability distribution shifts towards right as the mean value of work scales linearly with n. Fine 
structure in probability distributions get smeared out progressively and distribution tends towards a 
Gaussian. In the inset of figure [13] (a) we have shown the Gaussian fit for the obtained distribution 
for 10 cycles. From this fit we obtain (Wi 0p ) = 1.16 and variance a 2 = (W 2 ) — (W) 2 = 0.37. From 

this we can obtain dissipation ratio Rdi SS = 27wvp — 1> i- e -> ^ ne variance equals ^(W) which is 
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the required condition to satisfy SSFT when observed distribution is Gaussian [24J, |48|, |49J. The va- 
lidity of SSFT for work is also obviously seen from figure [13] (b) where we have plotted the symmetry 
function (In pP^\ ) versus f3W for work evaluated over different cycles as indicated in the figure. As 
we increase the number of periods from 1 to 10 the slope of symmetry function approaches 1. Hence 
our data suggests that the SSFT is satisfied even for finite number of periods. The number of periods 
above which SSFT is valid depends sensitively on the parameters in the problem. As already noted 
heat fluctuations over a cycle are large compared to work fluctuations. The heat fluctuations get an 
additional contribution from the internal energy (eq. (HI)). Even for a linear problem of harmonic 
oscillator, heat distribution measured over a single period does not satisfy SSFT as opposed to work 
fluctuations. The contribution from internal energy is supposed to dominate at very large values 



of Q, making the distribution P(Q) exponential in the large Q limit [15j, [16J. However, it may be 
noted that the distribution of the change in internal energy does not change with number of periods. 
Heat being an extensive quantity in time, distribution changes as we change the number of periods 
as shown in figure ITHf a) where we have plotted P(Q np ) for various values of n. As anticipated, 
by increasing n distribution tends towards the a Gaussian (see for n=10 cycles). The Gaussian fit 
for the P(Qnp) (inset of figurefH] (a)) gives the value for the variance as 0.56, and mean as 1.74. 
Thus dissipation ratio is 0.99, which is closer to unity, satisfying SSFT. In principle, one should be 
able to observe exponential tails for the distribution P(Q) in the large Q limit 16[. However, our 
simulations will not be able to detect it due to lack of required precision. As mentioned earlier, 
in the limit n — > oo, conventional SSFT holds for heat distributions j26|. In figure [14Tb). we have 
plotted the symmetry functions (In T^z^y) as a function of f3Q. The slope of the symmetry function 
approaches unity as we increase n, thereby suggesting the validity of SSFT. 



IV. CONCLUSION 



In conclusion, we have studied the nature of energy fluctuations in a biharmonically driven bistable 
system. This system is driven simultaneously with two periodic input signals of frequencies u and 2u, 
having phase difference between them. The presence of additional periodic drive induces particle 
confinement or localization in a preferred potential well. The degree of confinement analyzed in terms 
of the averaged value of the particle position (x) depends on the system parameters. We have shown 
that in spite of confinement, SR signal when quantified via the averaged work per period exhibits 
enhanced response. This is in sharp contrast to the case when confinement is induced by static tilt, 
which degrades SR. Surprisingly, the average input energy over a period is not very sensitive to cf) 
even though variation of (f) leads to significant particle pumping. However, changes in <fr does affect 
qualitatively the nature of hysteresis loop and distributions/fluctuations of work and heat. We have 
analyzed the fluctuations in work done, heat dissipated, and internal energy over a large but finite 
number of periods. Our data suggests that the SSFT for work and heat hold in this system. 
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